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I.  INTRODUCTION 


The  self-similar  solution  for  the  expansion  of  a  gas  into  a  vacuum  is 
well  known.1  In  the  present  study  the  latter  similarity  solution  is 
generalized  to  include  cases  where  the  gas  is  bound  by  a  shell  with  negli¬ 
gible  tensile  strength.  The  shell  is  assumed  to  envelop  the  gas  during  the 
entire  expansion  process.  An  exact  solution  is  obtained  in  quadrature 
form.  Analytic  expressions  for  dependent  variables  are  presented  for 
limiting  cases. 

An  approximate  local  similarity  solution  is  also  presented  which 
includes  the  effect  of  a  highly  compressible  porous  medium  external  to  the 
shell.  The  latter  solution  is  compared  with  the  solution  by  Zeldovich  and 

p 

Raizer  of  a  point  explosion  in  a  porous  medium. 
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I I .  THEORY 


We  consider  the  planar  (o  =  0),  cylindrical  (c  =  1)  and  spherical  (o  =•  2) 
case  of  a  gas  bound  by  a  shell-like  mass  which  is  surrounded  by  a  vacuum 
(Fig.  1).  The  expansion  is  assumed  to  start  at  time  t  =  0.  The  initial 
gas  radius  is  denoted  r0  and  the  corresponding  shell  location  is  denoted  Rq. 
The  initial  gas  conditions  at  r  =  0  are  also  denoted  by  subscript  zero.  We 
assume  that  the  entropy  of  the  gas  is  everywhere  the  same.  An  exact  simi¬ 
larity  solution  is  first  obtained.  An  approximate  local  similarity  solu¬ 
tion  which  includes  the  effect  of  a  porous  medium  external  to  the  shell  is 
then  indicated. 

A.  EXACT  SIMILARITY  SOLUTION 

The  equations  of  motion  for  the  gas  are 

|£  +  l£Y  £fiv  _  0  CONTINUITY  (la) 

3t  3r  r 

11  .  V  lx  .  .  i£  momentum  (1b) 

ot  3  r  y  3r 

p0 

where  p/pY  =  Pq/Pq*  Symbols  are  defined  in  Appendix  A.  The  location  of 
the  shell  at  each  instant  is  denoted  R(t).  We  introduce  the  similarity 
variable 

u  =  r/R  (2) 

so  that  the  gas  is  confined  to  the  region  -1  <  n  <  1  for  the  planar  case 
and  to  the  region  0  <  n  *  1  for  a  =  1,  2.  The  dependent  variables  are 
expressed  in  the  form 

„  p  -(o+1) 

ii 1  {k]  f(n) 

v  =  Rn  (3b) 
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(a)  PLANAR  (o=0) 


(b)  CYLINDRICAL  (o=1) 


(c)  SPHERICAL  (o=2) 


Fig.  1.  Initial  Conditions  for  Gas  Bound  by  Shell 
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which  satisfies  Eq.  (la).  Here,  R  =  dR/dt.  Substitution  into  Eq. 
yields 


-R 

*0<R/E0> 


R _ 

-( Y- 1 ) ( o+1 ) 


fY~2  df/dn  . 

n  ’  y-1 


where  C  is  a  constant.  Integration  of  Eq.  (4)  yields 


(1b) 


(4) 


£_  . 
p0 


5  1/(y-1) 
( 1-Cn  ) 


-(y-D(o+1) 


1/2 


(5a) 


(5b) 


where  R  is  the  speed  of  the  shell  after  long  times  and  equals 
00 

R  0  r  \/2 

—  =  —  f— 1 
aQ  y-1  '•o+V 

It  is  assumed,  in  Eq.  (5b),  that  R  =  0  at  t  =  0. 


(5c) 


The  constant  C  is  evaluated  by  considering  the  boundary  conditions  at 
n  =  1.  The  total  mass  of  the  shell  is  denoted  M„.  The  shell  mass,  per 
unit  area,  during  the  expansion  is 


SHELL  HASS  _  ^S_ 
AREA  =  {Ro 


where  6  =  2,  2*  and  4n  for  o  =  0,  1  and  2.  Note  that  gas  volume  is  given 
by  6R0+1/(o+1).  The  shell  is  accelerated  due  to  the  gas  pressure  at  n  =  1. 
The  appropriate  form  for  Newton's  law,  at  n  =  1,  is 

6R°p  =  M  R  (6) 

s  s 


where  ps  is  the  gas  pressure  at  n  =  1.  The  value  of  C,  which  is  consistent 
with  Eqs.  (5)  and  (6),  is 


(7) 


The  quantities  pQ  and  p0  can  be  related  to  the  net  gas  mass  Mg  and  the  net 


initial  gas  energy  E  via  the  expressions 

© 


( o+l )M  1  -  1 /( Y- 1 )  „ 

MI  =  - =  (<*♦!)  J  ( 1-Cn  )  - °dn 

*<>oC  0 


(8a) 


(o+1)(y-1)  E 


1 


El  = 


,  -  ■  o  y/(y~D  _ 

- — ^  =  (o+1)  J  (1-Cn  )  n  dn 

o 


(8b) 


Note  that  MI  equals  the  ratio  of  the  average  gas  density  at  t  =  0  to  pQ. 
Similarly,  El  equals  the  average  gas  energy  per  unit  volume  at  t  =  0  to  the 
initial  energy  per  unit  volume  at  the  origin,  Pq/(y-1).  Equation  (7)  can 
then  be  expressed 


Mg  2y  MI  C _ 

Ms  "  Y"1  0+1  (1-C)y/(t'1) 


which  relates  Mg/Mg  to  C  for  a  given  gas  and  geometry.  Numerical  values  of 

MI,  El  and  M_/M_  are  listed  in  Table  1  as  functions  of  y,  o,  C.  Hence, 

©  3 

Eqs.  (3)  and  (5)  can  be  evaluated  if  Yi  o,  and  Mg/Mg  are  specified.  It  is 

seen  that  0  <  C  <  1  for  0  <  M  /M  <  ®.  The  case  C  =  1  corresponds  to  an 

g  s 

unbound  gas. 

Analytic  expressions  for  MI  and  El  are,  respectively, 


MI 


o+ 1  C  o+ 1  2-y 

o+3  y-1  +  o+5  2(y-1)2 


C2  +  0  (C3) 


(9a) 
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Table  1.  Variation  of  Mass  Ratio  M_/Ms,  Normalized  Mass 

Integral  MI  and  Normalized  Energy  Integral  El  with 
11/9  <  y  <  5/3,  o  =  0,1,2  and  0.00  <  C  <  1.00 


SIGMA =0.0 

SIGMAsl .0 

SICMA=2 

.0 

GAMMA 

C 

MG/MS 

MI 

El 

MG/MS 

MI 

El 

MG/MS 

MI 

El 
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0.00 

0 . 000E+00 

1.0000 

1.0000 

0 . 000E*00 

1.0000 

1.0000 

0.000E+00 

0  ******** 

1  •  w w 

X  •  POIW 

0.02 

0 . 104E+00 

0.9900 
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0.618E-01 
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0.04 
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0.712E-01 
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0.06 

0 . 340E+00 
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0.40 
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0.6948 
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0.6933 

0.4371 
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0.4669 
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0.70 

0.490E+02 
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0.4911 

0.3669 
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MG/MS 

MI 
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0.2553 

0.1798 

0.90 

0.103E*0S 

0.6174 

0.4627 

0 . 316E+04 

0.3174 

0.2469 

0.143E*04 

0.2164 

0.1609 

1.00 

INFINITY 

0.4909 

0.4296 

INFINITY 

0.2857 

0.2222 

INFINITY 

0.1841 

0.1289 

GAMMA 

C 

MG/MS 

MI 

El 

MG/MS 

MI 

El 

MG/MS 

MI 

El 

1.286 

0.00 

0.000E+00 

1.0000 

1.0000 

0.000E+00 

1.0000 

x  •  trow 

0 . 000E+00 

X  •  WWW 

1.0000 

0.02 

0.193E+00 

0.9770 

0.9706 

0.952E-01 

0.9666 

0.9560 

0.630E-01 

0.9587 

0.9473 

0.04 

0.413E+00 

0.9647 

0.9425 

0 . 202E+00 

0.9323 

0.9141 

0 . 133E+00 

0.9189 

0.8973 

0.06 

0 .666E»00 

0.9331 

0.9155 

0 . 321E+00 

0.9001 

0.8741 

0.209E+00 

0.8806 

0.8497 

0.08 

0 . 956E+00 

0.9121 

0.8896 

0.465E+00 

0.8690 

0.8360 

0.29BE+00 

0.8436 

0.8046 

0.10 

0 . 129E*01 

0.8918 

0.8648 

0.607E+00 

0.8390 

0.7997 

0 . 389E+00 

0.8080 

0.7616 

0.20 

0 . 393E+01 

0.7992 

0.7659 

0 . 173E+01 

0.7040 

0.6426 

0. 106E+01 

0.6493 

0.6785 

0.30 

0 . 969E+01 

0.7206 

0.6686 

0 . 398E+01 

0.6919 

0.5208 

0.233E+01 

0.6197 

0.4400 

0.40 

0 . 236E+02 

0.6641 

0.6988 

0.896E+01 

0.4998 

0.4272 

0 . 496E*01 

0.4163 

0.3366 

0.60 

0.609E+02 

0.6983 

0.6429 

0.216E+02 

0.4248 

0.3666 

0. 113E+02 

0.3325 

0.2604 

0.60 

0 . 184E+03 

0.6617 

0.4981 

0.608E+02 

0.3644 

0.3011 

0.298E+02 

0.2677 

0.2048 

0.70 

0.728E+03 

0.6127 

0.4619 

0.224E+03 

0.3160 

0.2694 

0 . 103E+03 

0.2178 

0.1645 

0.80 

0.483E+04 

0.4801 

0.4322 

0. 140E+04 

0.2776 

0.2272 

0 . 603E+03 

0.1798 

0.1360 

0.90 

0. 118E+06 

0.4627 

0.4076 

0. 316E+06 

0.2469 

0.2020 

0.129E+06 

0.1609 

0.1132 

1.00 

INFINITY 

0.4296 

0.3866 

INFINITY 

0.2222 

0.1818 

INFINITY 

0.1288 

0.0966 

GAMMA 

C 

MG/MS 

MI 

El 

MG/MS 

MI 

El 

MG/MS 

MI 

El 

1.222 

0.00 

0.000E+00 

X  tlWW 

1.0000 

0.000E+00 

«  ******** 

1  .019190 

1.0000 

0.000E+00 

1  **^**^ 

X  •  ******** 

0  HUtUM 1 

x  •  ww 

0.02 

0 . 239E+00 

0.9706 

0.9643 

0.118E+00 

0.9660 

0.9466 

0.776E-01 

0.9473 

0.9361 

0.04 

0.619E*00 

0.9426 

0.9305 

0 . 262E+00 

0.9141 

0.8964 

0 . 165E+00 

0.8973 

0 • 8762 

0.06 

0 . 849E+00 

0.9166 

0.8986 

0.406E+00 

0.8741 

0.8491 

0.263E+00 

0.8497 

0.8201 

0.08 

0. 124E*01 

0.8896 

0.8681 

0 . 682E+00 

0.8360 

0.8046 

0.373E+00 

0.8046 

0.7676 

0.10 

0.170E+01 

0.8648 

0.8394 

0.785E+00 

0.7996 

0.7628 

0.499E«00 

0.7616 

0.7184 

0.20 

0.668E+01 

0.7669 

0.7173 

0.241E+01 

0.6426 

0.5888 

0 . 145E+01 

0.6786 

0.6174 

0.30 

0 . 167E+02 

0.6686 

0.6246 

0 .611E+01 

0.6208 

0.4623 

0 .344E+01 

0.4400 

0 • 3768 

0.40 

0 . 438E+02 

0.6987 

0.6539 

0.166E+02 

0.4271 

0.3707 

0.820E+01 

0.3366 

0.2773 

0.60 

0 . 136E+03 

0.6429 

0.4995 

0.443E+02 

0.3566 

0.3043 

0.216E*02 

0.2604 

0.2093 

0.60 

0 . 608E+03 

0.4981 

0.4671 

0.163E+03 

0.3010 

0.2667 

0.696E*02 

0.2048 

0.1623 

0.70 

0 . 267E*04 

0.4618 

0.4236 

0.761E+03 

0.2694 

0.2197 

0.317E+03 

0.1646 

0.1296 

0.80 

0 . 266E+06 

0.4322 

0.3962 

0.699E«04 

0.2272 

0.1923 

0.277E*04 

0.1360 

0.1061 

0.90 

0.128E+07 

0.4076 

0.3736 

0.317E+06 

0.2020 

0.1709 

0.118E+06 

0.1132 

0.0889 

1.00 

INFINITY 

0.3866 

0.3543 

INFINITY 

0.1818 

0.1538 

INFINITY 

0.0966 

0.0769 
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=  |1  -  (1  -  C)T/(T-1>I  (o=1)  (9b) 

EI.I-^-^tC.2^  - *— 5  C2  *  0(C3)  (10a) 

0.3  ,-1  o.5  2(y-1)2 

■  ^j}r!)c  I1  -  (i-c)(2t'1)/(t~1)1  (o=d  nob) 

Expressions  for  the  case  o  =  0,2  are  given  in  Appendix  B.  Thus,  H-/Ha  and 

g  s 

C  are  further  related  by 


JL  .  2-yC 
M  "  (y-1 )(o+1) 

5 


1 


(t 


o-*- 1 1  C 
o+3JY-1 


(11a) 


y/(y-1) 


-1 


(0=1)  (11b) 


The  shell  trajectory  can  be  expressed  in  the  form 


Rmt  R/Rq  d(R/R0) 

R0  =  J1  [1  -  (R/Ror<’,-')<0',)1,/2  ('2a) 

=  [(R/Rq)2-1 J1/2  (o=2,  r=5/3)  (12b) 


2((R/Rq)-  1]1/2 
[(y-1)(o+1)]1/2 


h  +  °  (f-  -  01 

n0 


(12c) 


=  (R/RQ)  [l  ♦  0  (R0/R)(y'1){o+1)] 


(12d) 


Corresponding  values  of  R^t/Rg  and  R/Rq  are  given  in  Table  2  for  various 
values  of  y  and  o.  The  acceleration  of  the  shell  decreases  with  time,  and 
the  shell  ultimately  moves  with  the  constant  velocity  given  by  Eq.  (5c). 


12 


Table  2.  Shell  Ordinate  R  =  R/RQ  Versus  Time  TAU  =  R^t/R.. 
for  11/9  s  y  -  5/3  and  o  s  0,1,2  “ 


TAU 


SIGMA s0.  SIGMA=1.  SIGMA=2. 


GAMMA 

GAMMA 

GAMMA 

GAMMA 

GAMMA 

GAMMA 

GAMMA 

GAMMA 

GAMMA 

GAMMA 

GAMMA 

GAMMA 

R 

6/3 

7/6 

9/7 

11/9 

5/3 

7/5 

9/7 

11/9 

6/3 

7/6 

9/7 

11/9 

1.0 

A  IMUUI 
*>  .  VXfW 

0.0000 

0.0000 

A  0iUUI 
v  •  WW 

A  AAA A 

V  •  WW 

A  AAAA 

V  .  Dvw 

A 

mf  •  vvw 

0.0000 

A  AAAA 
mJ  •  Vvw 

A  AAAA 
mf  •mfXfW 

A  AAAA 
V  •  wW 

A  AAAA 

Wf  .  ww 

1.6 

1.8442 

2.3674 

2.7774 

3.1419 

1.3365 

1.6919 

1.9850 

2.2402 

1.1180 

1.4019 

1.6380 

1.8443 

2.0 

2.7494 

3.4869 

4.0927 

4.6206 

2.0322 

2.6326 

2.9608 

3.3172 

1.7320 

2.1236 

2.4560 

2.7496 

2.6 

3.6264 

4.4390 

6.1968 

6.8562 

2.6487 

3.2583 

3.7743 

4.2291 

2.2912 

2.7690 

3.1646 

3.6266 

3.0 

4.2402 

6.3070 

6.1954 

6.9723 

3.2304 

3.9304 

4.6300 

6.0612 

2.8284 

3.3560 

3.8219 

4.2406 

3.6 

4.9184 

6.1231 

7.1313 

8.0149 

3.7924 

4.6707 

6.2448 

6.8448 

3.3540 

3.9310 

4.4493 

4.9187 

4.0 

6.6720 

6.9038 

8.0233 

9.0065 

4.3419 

6.1899 

6.9318 

6.6954 

3.8729 

4.4918 

6.0668 

6.6723 

4.6 

6.2078 

7.6683 

8.8828 

9.9604 

4.8828 

6.7938 

6 . 6986 

7.3216 

4.3874 

6.0428 

6.6601 

6.2081 

6.0 

6.8299 

8.3927 

9.7172 

10.8849 

6.4174 

6.3863 

7.2600 

8.0291 

4.8989 

6.6866 

6.2326 

6.8303 

6.6 

7.4413 

9.1109 

10.6313 

11.7857 

5.9474 

6.9697 

7.8891 

8.7217 

6.4083 

6.1247 

6.8068 

7.4417 

6.0 

8.0440 

9.8168 

11.3287 

12.6669 

6.4736 

7.6469 

8.6181 

9.4019 

6.9160 

6.6686 

7.3741 

8.0444 

6.6 

8.6394 

10.6096 

12.1119 

13.6316 

6.9968 

8.1160 

9.1388 

10.0719 

6.4226 

7.1888 

7.9360 

8.6398 

7.0 

9.2287 

11.1937 

12.8831 

14.3819 

7.6177 

8.6811 

8.7624 

10.7331 

6.9281 

7.7162 

8.4932 

9.2291 

7.6 

9.8127 

11.8696 

13.6437 

16.2199 

8.0366 

9.2419 

10.3698 

11.3866 

7.4330 

8.2412 

9.0464 

9.8132 

8.0 

10.3921 

12.6383 

14.3960 

16.0469 

8.6637 

9.7990 

10.9619 

12.0336 

7.9372 

8.7641 

9.6962 

10.3926 

8.6 

10.9676 

13.2006 

16.1381 

16.8642 

9.0694 

10.3628 

11.6693 

12.6746 

8.4409 

9.2863 

10.1430 

10.9681 

9.0 

11.6394 

13.8671 

16.8739 

17.6728 

9.6839 

10.9037 

12.1626 

13.3102 

8.9442 

9.8049 

10.6872 

11.6399 

9.6 

12.1081 

14.6086 

16.6029 

18.4735 

10.0973 

11.4521 

12.7420 

13.9413 

9.4471 

10.3232 

11.2291 

12.1086 

10.0 

12.6739 

16.1662 

17.3269 

19.2669 

10.6097 

11.9982 

13.3281 

14.6680 

9.9498 

10.8403 

11.7689 

12.6745 

Flow  profiles,  the  net  impulse  on  the  plane  of  symmetry  at  r  =  0,  and 
flow  conditions  at  r^>  rQ  are  now  noted. 

Flow  Profiles:  Equation  (5)  can  be  expressed 


r  0+1  6Rn+1 

(r)  0 


_  0  .  11  -  cn  ) 

(o+1)Mg  p  ‘  MI 


1 / ( Y~ 1 ) 


(R_jTto*1>  6R0+1 

V 


_ .  _E_  .  O-Cn  ) 

(o+1 )E  y-1  “  El 

© 


y/(y-1) 


(13a) 


(13b) 


where  p  and  p  have  been  normalized  in  terms  of  Mg  and  Eg,  respectively. 
Equation  (13a)  is  plotted  in  Fig.  2.  Density  and  pressure  are  uniform, 
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t 


with  ti,  for  the  case  C  =  0.  In  the  latter  case  the  shell  mass  is  large  and 
moves  slowly  (R  ♦  0) ,  so  that  gas  density  and  pressure  are  equilibrated 

<D 

in  the  region  0  <  n  s  1.  The  density  and  pressure  decrease  at  n  =  1 ,  as  C 

increases,  and  ultimately  become  zero  when  C  =  1 .  As  previously  noted,  the 

case  C  =  1  corresponds  to  the  expansion  into  vacuum  of  an  unbound  gas. 

Impulse  at  r  =  0:  The  net  impulse  on  the  plane  of  symmetry  at  r  =  0  is  de¬ 
noted  Iq  and  is  found  from 

o>  A 

In  =  J  dt  J  pdA  (14) 

u  0  0 

2 

where  A  =  1,2R,  and  nR  for  o  =  0,1,  and  2,  respectively.  The  result  may 
be  expressed 


where 


R 


2B 

(y-D(o-H) 


B  =  1/2 


■i  ' 


(1 


9  y/(y-1) 

Cn  )  dn 


=  4C  2y-1 


[1  -  (l-C)(2*-1)/<*-1)] 


(15a) 

o=0  (15b) 

o=1  (15c) 

o=2  (15d) 


and  (El)  n  denotes  El  evaluated  with  o  =  0.  Normalization  in  terms  of  gas 
o=u 

mass  and  initial  gas  energy  yields 


For  C  ♦  0 


<2EgMe> 


1/2 


’  '2yC  MI  El  ' 


M  1/2  I0 

‘"s'  (2EM)’/2 

g  g 


*£,£,][  ^  o  =  0,1,2 


(16a) 


(16b) 
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Equation  ( 16a)  is  listed  in  Table  3.  The  net  impulse  lQ/(2EgMg)^ 
increases  with  decrease  in  C  because  the  shell  mass  impedes  the  expansion 
and  thus  the  pressure  at  r  =  0  remains  relatively  high  for  a  longer  period 
of  time.  In  the  limit  C  +  0,  the  initial  gas  internal  energy  is 
transferred  to  shell  kinetic  energy. 


Flow  at  r  >>  rQ:  In  the  limit  r/rQ  >>  1,  R+  R^  and  Eqs.  (3)  and  (5)  become 

v  =  Rb  n  (17a) 

R  =  R  t  (17b) 

CD 

Let  subscript  1  denote  conditions  at  r1  >>  rQ.  The  variation  of  gas 
dynamic  pressure  with  time  at  this  station  is 


(R1/Rn)0+1p.v  2  t  o+3  ,  t.  - 

-LJL  » (r)  1  -  C  C^)2] 

p0r;  fc 


1/( Y-1 ) 


(18) 


where  t^  =  r^/R^  and  t^/t  <  1.  The  maximum  dynamic  pressure  is  denoted 

(p,v,2)  and  occurs  at  time  t_.  These  can  be  expressed  in  the  form 
l  i  m  “• 


and 


m 


g 


4yC 

(y-I)EI 


t.  ?  1/(y-1) 

dr)  1 

cm 


(19a) 


—  =  G  G  <  1  (19b) 

m 

=1  G  >  1  (19c) 
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Table  3 


Net  Impulse  at  Plane  of  Symmetry  IM  =  Iq/(2M  E  )1/2 
for  11/9  <  y  <  5/3,  o  =  0,1,2  and  0.02  <  C  I  1.00. 


IM 


GAMMA 

C 

SIGMAs0.0 

SIGMAsl .0 

SXGMAs2 

1.667 

0.00 

00 

00 

00 

0.02 

1.602 

1.428 

1.368 

0.04 

1.148 

1.014 

0.967 

0.06 

0.950 

0.831 

0.789 

0.08 

0.834 

0.722 

0.683 

0.10 

0.766 

0.649 

0.610 

0.20 

0.571 

0.469 

0.431 

0.30 

0.499 

0.394 

0.353 

0.40 

0.461 

0.352 

0.309 

0.50 

0.440 

0.326 

0.280 

0.60 

0.428 

0.310 

0.261 

0.70 

0.421 

0.301 

0.249 

0.80 

0.417 

0.295 

0.242 

0.90 

0.416 

0.293 

0.238 

1.00 

0.416 

0.292 

0.238 

1.400 

0.00 

00 

00 

00 

0.02 

1.363 

1.211 

1.169 

0.04 

0.983 

0.863 

0.820 

0.06 

0.819 

0.709 

0.671 

0.08 

0.723 

0.619 

0.582 

0.10 

0.660 

0.558 

0.522 

0.20 

0.513 

0.412 

0.374 

0.30 

0.459 

0.363 

0.312 

0.40 

0.434 

0.323 

0.278 

0.60 

0.422 

0.306 

0.258 

0.60 

0.416 

0.297 

0.246 

0.70 

0.413 

0.292 

0.240 

0.80 

0.412 

0.291 

0.237 

0.90 

0.412 

0.290 

0.236 

1.00 

0.412 

0.290 

0.236 

1.286 

0.00 

00 

00 

00 

0.02 

1.210 

1.072 

1.024 

0.04 

0.879 

0.766 

0.726 

0.06 

0.736 

0.632 

0.595 

0.08 

0.654 

0.654 

0.518 

0.10 

0.600 

0.601 

0.465 

0.20 

0.479 

0.377 

0.338 

0.30 

0.438 

0.330 

0.287 

0.40 

0.421 

0.307 

0.261 

0.60 

0.414 

0.296 

0.247 

0.60 

0.410 

0.291 

0.240 

0.70 

0.409 

0.289 

0.236 

0.80 

0.409 

0.289 

0.235 

0.90 

0.409 

0.289 

0.235 

1.00 

0.409 

0.289 

0.236 

1.222 

0.00 

00 

00 

00 

0.02 

1.102 

0.973 

0.928 

0.04 

0.806 

0.698 

0.660 

0.06 

0.679 

0.578 

0.642 

0.08 

0.606 

0.608 

0.473 

0.10 

0.669 

0.461 

0.426 

0.20 

0.458 

0.354 

0.314 

0.30 

0.426 

0.316 

0.271 

0.40 

0.414 

0.299 

0.260 

0.50 

0.409 

0.291 

0.241 

0.60 

0.408 

0.289 

0.236 

0.70 

0.407 

0.288 

0.235 

0.80 

0.407 

0.288 

0.236 

0.90 

0.407 

0.288 

0.234 

1.00 

0.407 

0.288 

0.234 
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where 


_ !_  [  (o-*-3)(y-1)1/^ 

•  1/2  l(o+3)(y-1)  ♦  2‘ 


09d) 


In  the  limit  C  ♦  0, 


t1/tm  = 


(20a) 


6R 


o+l 


(»,v,2) 


m 


g 


Hr 

M 

s 


2(o+1)  r*  [1  +  0(C)] 


(20b) 


p 

Numerical  results  for  (p,v,  )  and  t_  are  given  in  Table  4.  Values  of 

l  l  m  m 

t^tj,  in  the  range  1 1/tm  <  1  indicate  a  local  maximum  in  dynamic 
pressure.  The  case  t^t^  =  1  corresponds  to  a  flow  wherein  the  maximum 
dynamic  pressure  occurs  directly  behind  the  shell.  The  maximum  dynamic 
pressure  is  seen  to  decrease  with  decrease  in  C  (i.e.,  increase  in  shell 
mass).  The  latter  trend  occurs  because,  at  a  fixed  location  r^,  the  gas 
density  is  relatively  unaffected  by  the  shell  whereas  the  flow  velocities 
are  decreased  due  to  shell  inertia. 


The  impulse  per  unit  area  at  r1 ,  due  to  gas  flux,  is  denoted  Ig  and  is 
found  from 


Table  4. 


Maximum  Dynamic  Pressure  MDP  =  6R?+\p.v  ^ )  /E 
Corresponding  Time  XM  =  t1/tm  for  r,  >>  r0,® 
11/9  Z  y  <  5/3,  a  =  0,1,2  and  0.02  <  C  <  1.00 


g 


SIGMA=0.0 


SIGMA=1 .0  SIGMA=2 .0 


GAMMA  C  XM  MDP 


XM  MDP 


XM  MDP 


1.667  6.00  1.0000 

6.02  1.0000 
0.04  1.0000 

0.06  1.0000 
0.08  1.0000 
0.10  1.0000 
0.20  1.0000 
0.30  1.0000 

0.40  1.0000 

0.S0  1.0000 

0.60  0.9129 

0.70  0.8452 

0.80  0.7906 

0.90  0.7454 

1.00  0.7071 


1.400 


0.00 

0.02 

0.04 

0.06 

0.08 

0.10 

0.20 

0.30 

0.40 

0.60 

0.60 

0.70 

0.80 

6.90 

1.00 


1.0000 

1.0000 

1.0000 

1.0000 

X  •  1/vUU 

1*^0 
•  vuvv 


0.9682 

0.8680 

0.7906 

0.7319 

0.6847 

0.6455 

0.6124 


1.286  0.00  1.0000 
0.02  1.0000 
0.04  1.0000 

0.06  1.0000 
0.08  1.0000 
0.10  1.0000 
0.20  1.0000 
0.30  1.0000 

0.40  0.8660 

0.60  0.7746 

0.60  0.7071 

0.70  0.6646 

0.80  0.6124 

0.90  0.6773 

1.00  0.6477 


1.222  0.00  1.0000 
0.02  1.0000 
0.04  1.0000 

0.06  1.0000 
0.08  1.0000 
0.10  1.0000 
0.20  1.0000 
0.80  0.9128 

0.40  0.7906 

0.60  0.7071 

0.60  0.6466 

0.70  0.5976 

0.80  0.6590 

0.90  0.6270 

1.00  0.6000 


0.0000 

0.1973 

0.3890 

0.6748 

0.7644 

0.9277 

1.6876 

2.2453 

2.6688 

2.6338 

2.5833 

2.6591 

2.6494 

2.6468 

2.6466 

0.0000 

0.2726 

0.6297 

0.7712 

0.9969 

1.2064 

2.0056 

2.3896 

2.3998 

2.3466 

2.3234 

2.3146 

2.3120 

2.3116 

2.3116 

0.0000 

0.3466 
0.6623 
0.9500 
1.2090 
1.4396 
2 . 1809 
2.3178 
2.2413 
2.2109 
2 . 1998 

2.1966 
2.1968 

2.1967 

2.1967 

0.0000 

0.4167 

0.7871 

1.1122 

1.3931 

1.6313 

2.2471 

2.2026 

2.1611 

2.1333 

2.1280 

2.1267 

2.1266 

2.1266 

2.1265 


•  oooo 

•  oooo 
.  oooo 

1 . 0000 
x  •  row 

1.0000 

I.... 

.row 

ICkOMOk 

.roro 

1.... 

•  roro 

1  ******** 
•roro 

0.9769 

0.9036 

0.8462 

0.7968 

0.7659 

1.0000 

1  ******** 

•  roro 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.9428 

0.8607 

0.7968 

0.7464 

0.7027 

0.6667 

1.0000 
1.0000 
1.0000 
x  •  roro 
1.0000 
1.0000 
1.0000 
1.0000 
0.9634 
0.8628 
0.7785 
0.7207 
0.6742 
0.6366 
0.6030 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.8770 

0.7844 

0.7161 

0.6630 

0.6202 

0.6847 

0.6647 


0.0000 

0.1989 

0.3966 

0.6897 

0.7809 

0.9689 

1.8480 

2.5873 

3.1266 

3 . 3936 

3.3418 

3.2646 

3.2180 

3.2076 

3.2066 

0.0000 

0.2767 

0.6424 

0.7996 

1.0463 

1.2822 

2.2766 

2.9088 

3.1246 

2.9962 

2.9099 

2.8756 

2.8649 

2.8629 

2.8628 

0.0000 

0.3609 
0.6828 
0.9960 
1.2866 
1.5568 
2.5654 
2.9763 
2.8636 
2.7619 
2 . 7087 
2.6947 
2.6915 
2.6911 
2.6911 

0.0000 

0.4246 

0.8170 

1.1769 

1.6031 

1.7962 

2.7376 

2.8675 

2.6847 

2.6166 

2.6944 

2.6887 

2.6877 

2.6877 

2.6877 


1.0000 

1  ******** 
*  roro 

1  ******** 
.  row 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.9449 

0.8839 

0.8333 

0.7906 

1.0000 
x .  row 
1.0000 
1.0000 
1.0000 

1  ******** 

.  roro 

1.0000 

1  ******** 

.  roro 

1.0000 
x  .roro 
0.9129 
0.8462 
0.7906 
0.7454 
0.7071 

1.0000 

1  ******** 

.  roro 

1.0000 

1.0000 
,  mt*^** 

x  .roro 
1.0000 
1.0000 
1.0000 
1.0000 
0.9129 
0.8333 
0.7716 
0.7217 
0.6804 
0.6466 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.9449 

0.8461 

0.7716 

0.7143 

0.6681 

0.6299 

0.5976 


XI  ******** 

v .  row 

0.2000 

0.3997 

0.5990 

0.7976 

0.9962 

1.9669 

2.8366 

3.6641 

4.0443 

4.1672 

3.9723 

3.8823 

3.8663 

3.8526 

0  ******** 

.  roro 

0.2777 

0.6603 

0.8172 

1.0778 

1.3314 

2.4686 

3.3133 

3.7600 

3.7217 

3.6160 

3.4296 

3.4011 

3.3966 

3.3963 


0.0000 
0.3541 
0.6966 
1.0236 
1.3369 
1.6346 
2.8496 
3.6221 
3.6790 
3.3214 
3.2121 
3.1747 
3 . 1666 
3.1644 
3 . 1644 

0.0000 

0.4292 

0.8369 

1.2186 

1.6767 

1.9061 

3.1167 

3.6277 

3.2726 

3.1025 

3.0436 

3.0274 

3.0244 

3.0242 

3.0242 
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Integration  yields 


“06! 


gi  [1  -  (i-o’r/(l-1)| 


*  5  m)o-.! 


o=0  (21a) 
o=1  (21b) 


(  -Y-J — 

l(2Y-1)C 


[,  .  (,.c)(2,-1)/(r-D| 


-  (1-C)T/<,'')1 


0=2 


(21c) 


where  (MI)o_2  indicates  that  MI  is  evaluated  with  o  =  2.  The  impulse  per 
unit  area  at  due  to  the  shell  is  Ig  =  ^/(fiR^  ),  which  can  be 

expressed 


<VR0>°*'  l3 


P0R-  fc1 


(22) 


The  expressions  for  Ig  and  Is  can  be  put  in  the  form 


wherein  Ig  and  6R1°  Ig  are  referenced  to  the  net  impulse  associated 
with  a  free  gas  expansion.  In  the  limit  C  ♦  0,  Eqs.  (23a)  and  (23b)  become 


1/2  6R°  I 


(2M  E  ) 
g  g 


i  M 

£ _  0+2  _£ 

1/2  ‘  o+2  M 

s 


(2Ha) 


1/2 


fiR1°  1 s 

(2M  E  )1/2 

g  g 


(2^b) 


#  p 

In  the  limit  C  ♦  0,  the  shell  kinetic  energy  Mjl^/2  equals  the  original 
gas  internal  energy  E_  and  the  quantity  6R?I  corresponds  to  the  net 

o  IS 

momentum  of  a  mass  Mg  with  kinetic  energy  Eg. 

Numerical  results  for  I_,  I_  and  I.  s  I  ♦  I  are  given  in  Table  5. 

n  ®1/P3  n  ®  S  i /o 

The  quantities  6R,  I./(2M  E  ),/£  and  fiR,  I  /(2M  E  ) '  c  increase 

O  It  -I  /o6  g  '  ®  g  g 

while  6R-  I  /  (2M  E  )  decreases  as  C  decreases.  This  behavior  is  again 

l  g  g  g 

due  to  the  reduction  in  gas  expansion  rate  as  C  decreases. 

B.  APPROXIMATE  LOCAL  SIMILARITY  SOLUTION 

The  problem  is  now  generalized  to  include  a  highly  compressible  porous 
medium  external  to  the  shell.  The  configuration  is  illustrated  in  Fig.  3. 
As  the  shell  expands,  the  porous  medium  is  compressed  and  adheres  to  the 
shell  in  the  form  of  a  thin  band  (i.e.,  an  inelastic  collision  is  assumed). 
The  subsequent  motion  is  deduced  herein  using  a  local  similarity 
approximation. 

Let  Mg  q  denote  the  mass  of  the  shell  at  t  =  0,  let  Mg  denote  the  mass 
of  the  shell  at  subsequent  times,  and  let  denote  the  density  of  the 
porous  medium.  The  shell  mass  at  any  instant  can  be  expressed 

M  M  . 

5s  =  -S*2  +  o(Y-1)  (25a) 
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Table  5.  Impulse  per  Unit  Area  at  r^  >>  r0  Due  to  Gas  IMG  =  6R°I  /(2M  E  )1/2, 
Due  to  Shell  IMS  =  6R1°Is/(2MgEg) 1/2  and  Due  to  Sum  IMT8=  IMG  + 

IMS  as  a  Function  of  11/9  S  y  <  5/3,  a  =  0,1,2  and  0.02  <  C  <  1.00. 


SIGMA=0.0 


SIGMA=1 .0 


GAMMA  C  IMG 


IMS  IMT  IMG 


IMS  IMT 


1.607  0.00  0.0000 

0.02  0 . 1578 

0.04  0.2228 

0.06  0.2723 

0.08  0.3136 

0.10  0.3S01 

0.20  0.4886 

0.30  0.6883 

0.40  0.6649 

0.60  0.7242 

0.60  0.7688 

0.70  0 . 8002 

0.80  0 . 8200 
0.90  0.8296 

1.00  0.8317 


oo  0 
3.0469  3.2048  0 
2.0738  2.2966  0 
1.6280  1.9003  0 
1.3640  1.6679  0 
1.1618  1.6118  0 
0.6642  1.1428  0 
0.4088  0.9971  0 
0.2671  0.9221  0 
0.1566  0.8797  0 
0.0866  0.8663  0 
0.0416  0.8417  0 
0.0149  0.8349  0 
0.0026  0.8322  0 
0.0000  0.8317  0 


0000 

00 

00 

.1494 

4.3381 

4.4876 

,2116 

2.9728 

3.1844 

,2697 

2.3501 

2 . 6097 

3003 

1.9685 

2.2689 

3364 

1.7012 

2.0376 

4788 

0.9949 

1.4737 

6886 

0.6477 

1.2362 

6792 

0.4258 

1.1050 

7650 

0.2700 

1.0250 

8169 

0.1579 

0.9748 

8646 

0.0797 

0.9442 

8970 

0.0303 

0.9273 

9142 

0.0056 

0.9199 

9183 

0.0000 

0.9183 

0.00 

0.0000 

oo 

0.02 

0.1861 

2 . 5403 

0.04 

0.2618 

1.7048 

0.06 

0.3189 

1.3188 

0.08 

0.3661 

1.0804 

0.10 

0.4069 

0.9124 

0.20 

0.6666 

0.4701 

0.30 

0.6663 

0.2637 

0.40 

0.7233 

0.1463 

0.60 

0.7689 

0.0746 

0.60 

0.7976 

0.0336 

0.70 

0.8136 

0.0122 

0.80 

0.8208 

0.0029 

0.90 

0.8229 

0.0003 

1.00 

0.8231 

0.0000 

00 

0.0000 

00 

2.7285 

0.1764 

3.6290 

1.9666 

0.2494 

2.4605 

1.6377 

0.3053 

1.9235 

1.4466 

0.3623 

1.6925 

1.3193 

0.3936 

1.3696 

1.0266 

0.6624 

0.7413 

0.9190 

0.6674 

0.4417 

0 . 8686 

0.7640 

0.2594 

0.8435 

0.8196 

0.1420 

0.8312 

0.8646 

0.0685 

0.8267 

0.8923 

0 .0265 

0.8237 

0.9061 

0.0068 

0.8232 

0.9107 

0.0006 

0.8231 

0.9111 

0.0000 

O0 

3 . 8064 
2.7098 
2.2288 
1.9448 
1.7631 
1.2937 
1.1091 
1.0140 
0.9616 
0.9330 
0.9188 
0.9129 
0.9113 
0.9111 


0.00 

0.0000 

00 

00 

0.02 

0.2103 

2.2100 

2.4204 

0.04 

0.2949 

1.4621 

1.7570 

0.06 

0.3578 

1.1146 

1.4724 

0.08 

0.4093 

0.8990 

1 . 3083 

0.10 

0.4632 

0.7471 

1.2003 

0.20 

0.6076 

0.3513 

0.9689 

0.30 

0.7007 

0.1761 

0.8768 

0.40 

0.7676 

0.0846 

0.8421 

0  ■  60 

0.7906 

0.0366 

0.8271 

0 . 60 

0.8076 

0.0133 

0.8209 

0.70 

0.8151 

0.0036 

0.8187 

0.80 

0.8176 

0.0006 

0.8181 

0.90 

0.8180 

0.0000 

0.8181 

1.00 

0.8180 

0.0000 

0.8180 

0.0000 

0.1996 

0.2816 

0.3439 

0.3969 

0.4412 

0.6106 

0.7260 

0.8040 

0.8656 

0.8856 

0.9001 

0.9054 

0.9065 

0.9066 


00 

3.1678 

2.1247 

1.6423 

1.3436 

1.1327 

0.6741 

0.3114 

0.1619 

0.0768 

0.0297 

0.0087 

0.0016 

0.0001 

0.0000 


00 

3.3674 
2.4063 
1 . 9862 
1.7395 
1.6739 
1 . 1846 
1.0363 
0.9669 
0.9313 
0.9163 
0.9089 
0.9069 
0.9066 
0.9066 


IMG 

0.0000 

0.1376 

0.1952 

0.2400 

0.2782 

0.3123 

0.4498 

0.6602 

0.6569 

0.7403 

0.8136 

0.8741 

0.9186 

0.9440 

0.9605 

0.0000 
0.1624 
0.2303 
0.2827 
0.3272 
0 . 3666 
0.5228 
0.6426 
0.7394 
0.8170 
0.8755 
0.9148 
0.9361 
0.9437 
0.9446 

0.0000 
0.1839 
0 . 2604 
0.3192 
0.3688 
0.4126 
0.5820 
0.7053 
0.7979 
0.8641 
0.9067 
0.9296 
0.9386 
0.9406 
0.9407 


SIGMA=2 .0 


IMS 

IMT 

oo 

00 

6.3346 

6.4721 

3.6709 

3.8661 

2.9144 

3.1644 

2.4619 

2.7301 

2.1284 

2.4407 

1.2763 

1.7251 

0.8634 

1.4136 

0.6789 

1.2348 

0.3803 

1.1206 

0.2316 

1.0451 

0.1221 

0 . 9962 

0.0486 

0.9672 

0.0096 

0.9536 

0.0000 

0.9505 

00 

oo 

4.4718 

4.6342 

3.0609 

3.2812 

2.4004 

2.6832 

2.0006 

2.3278 

1.7196 

2.0862 

0.9729 

1.4957 

0.6046 

1.2470 

0.3721 

1.1115 

0.2147 

1.0317 

0.1094 

0.9849 

0.0449 

0.9597 

0.0122 

0.9484 

0.0012 

0.9449 

0.0000 

0.9446 

oo 

00 

3.9114 

4.0963 

2.6466 

2.9059 

2.0626 

2.3818 

1.7025 

2.0713 

1.4486 

1.8610 

0.7717 

1.3536 

0.4428 

1.1481 

0.2461 

1.0430 

0.1226 

0.9867 

0.0615 

0.9582 

0.0162 

0.9457 

0.0030 

0.9416 

0.0001 

0.9407 

0.0000 

0 . 9407 
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where 


a  s 


»0  i  B0 

( o+1 )  H 


0+  1 


g 


Y  =  (R/R0) 


_J_  !e 

(MI )q  Pq 


o+l 


(25b) 


(25c) 


and  (MI)q  denotes  the  value  of  MI  at  t  =  0.  Conservation  of  momentum  can 
be  expressed  in  the  form 

(26) 


6R  p  =  d(M  R)/dt 
s  s 


Integration  yields 


(MSR) 

2M  E 
g  g 


(El) 


Y  M  p 

-  ;  /^dY 

0  1  Mg  P0 


(27) 


where  Eg  again  denotes  net  internal  gas  energy  at  t  =  0  and  where  (EI)0 
denotes  the  value  of  El  at  t  s  0.  In  order  to  integrate  Eq.  (27),  it  is 
necessary  to  express  Ps/Pg  as  a  function  of  Mg/Mg  or  Y.  The  assumption  of 
local  similarity  is  now  introduced.  Namely,  it  is  assumed  that  at  each 
instant  the  local  fluid  property  values  equal  those  for  a  self-similar  flow 
with  the  same  values  of  Ms/Mg  and  Y.  Thus,  recalling  Eqs.  (5a)  and  (8a), 
flow  properties  are  given  by 


■G-) 


1/y 


1  -  Cn2  (MI)q 
Y  MI 


v/R  =  n 


and  Ps/Pg  ca°  be  expressed  in  the  form 


£a 

p0 


4>Y"Y 

1  +  [ (o+l )/2] (M/M) 

5  S 


(28a) 
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where 


4> 


yC  MI 
Y-1 


a  M  (MI )n  Y 

I'^rl  I-mt) 

g 


=  i  [ (HI)0/MI 1Y 


(28b) 


For  o=1 


i>  =  i 


Y 

The  dependence  of  ip  on  M_/M_ 

o  6 

Ip  = 


1  ♦  (M  /M  )  -  [l  (M  /M)]1/y 
g  s  g  s  J 

,  for  o=0  and  2,  has  the  limiting  values 


1 


M  /M  =  0  (C  =  0) 
g  s 

M  /M  =  0  (C  =  1) 
^  © 


It  is  seen  that  t  is  an  order  1  quantity  and  has  a  weak  dependence  on 
Mg/Mg.  Thus,  the  major  dependence  of  Ps/Pq  on  Mg/Ms  is  explicitly 
displayed  in  Eq.  (28a).  Substitution  into  Eq.  (27)  yields 


(M  R)‘ 

2mV 

g  g 


(El) 


0  1 


Y  (M  /M 

I  “JLJL 


B  +  oY 


♦Y"y  dY 


(29) 


where 


B 


g 


-a 


Equation  (29)  can  be  integrated  to  obtain  MR  as  a  function  of  Y.  The 

9 

dependence  of  *  on  Mg/Mg  is  found  from  Eq.  (28b).  The  shell  trajectory  is 
found  from 
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(30) 


1 17 

(2h:/2t_  rR/R0  Ms  <2w  dR 

1  MgJ  R0  =  J1  Mg  MsR  R0 


Limiting  cases  are  now  noted.  In  the  limits  aY  <<  B,  oY  >>  B,  and 
B  =  0,  Eq.  (29)  becomes,  respectively, 


(E!)0  (MsR)‘ 


Y-1  2M  E 
g  g 


1 

B 


1 

a 


a 


/  ♦  Y'Y  |l  -§*♦  0(|^)2]  dY 

1  g 


/♦  *-T-’  £)  h  -  fe  *  0(|L)2]  dY 
i  g 


y  ,  K 

J  4  y"Y'  (jr)  dY  (B  =  0) 
1  g 


(31a) 


(31b) 


(31c) 


These  can  be  readily  integrated  analytically  if  a  mean  value  of  4  is 
used.  Equation  (31a)  treats  the  case  where  the  porous  medium  introduces 
only  small  departures  from  a  self-similar  flow.  Equation  (31b)  treats  a 
porous  medium  at  late  times.  Equation  (31c)  provides  a  particular  case 
which  is  valid  for  all  values  of  aY. 


In  the  limit  aY 
(31a)  yields 


0,  M  M  the  flow  is  self-similar. 
’  s  s,0’ 


(M..o/w.)‘ 


<EI)o  \omt  *  (0*,)/2 


Equation 


(32) 


which  agrees  with  Eq.  (5b).  The  role  of  0/Mg  is  explicitly  displayed  in 
Eq.  (32)  except  for  the  weak  dependence  of  4  on  MSj0/Mg. 

The  limit  aY  ♦  »  corresponds  to  the  case  of  a  porous  medium  and  large 
distances  from  the  origin.  In  this  limit  the  shell  trajectory,  velocity, 
momentum,  and  kinetic  energy  are,  respectively, 
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h-  (J 


>•*•1)  + 
4 


^  (MI)J  2E ^t2 
2~y  (EI)0  aM  R2^ 


(33a) 


M  1/2 


*  »  i  /  *-  -  V  *  *+  /  A  i 

(~§— )  R  =  [i~  Tpr r~  - ] 

(2EJ  l2-y  (EI)0  oyYJ 


(33b) 


(2M  E  ) 
g  g 


-  !*=■ 
*2- 


1 =1  1!1^0  >v2-Y1 

2-y  (EI)n 
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Vl  y-i  (mi)oT  1 
2Eg  =  2-y  (EI)0  yY-1 


(33d) 


It  is  seen  that  shell  radius  and  shell  momentum  increase  with  Y,  whereas 
shell  velocity  and  shell  kinetic  energy  decrease  with  Y.  The  decrease  in 
kinetic  energy  is  due  to  the  assumption  of  inelastic  collisions  between  the 
shell  and  the  porous  medium. 

It  is  also  seen  from  Eqs.  (32)  and  (33b)  that  when  the  expanding  gas 
is  bound  by  a  shell  and  a  porous  medium,  the  shell  velocity  increases  to 
some  maximum  value  and  then  decreases  during  the  expansion  process.  This 
behavior  can  be  illustrated  by  considering  the  limit  Mg/Mg  g  «  1,  a  = 

0(1),  In  this  limit  (EI)0  =  0  =  1,  B  =  Mg>0/Mg  and  Mg/Mg  =  (Mg>0/Mg) 

+  oY.  Equation  (29)  becomes 


(MR)4 

2ME 
g  g 


=  h _ — ]  +  [y2-y-i1 

H  l1  VY~ 1  2-y  11  1J 


The  first  term  in  the  right-hand  side  of  Eq.  (34)  corresponds  to  a  shell 
without  a  porous  medium  and  agrees  with  Eq.  (32)  in  the  limit  Mg/Ms  Q  •*  0. 
The  second  term  in  the  right-hand  side  of  Eq.  (34)  indicates  the  effect  of 
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a  porous  medium  and,  in  the  limit  aY  -*•  ®,  agrees  with  Eq.  (33b).  Equation 
(3*0  provides  a  solution  which  is  valid  for  all  values  of  Y. 

A  solution  for  a  strong  explosion  in  an  infinite  porous  medium  is  pre¬ 
sented  in  Ref.  2.  That  solution  assumes  that  the  decrease  in  kinetic 
energy  of  the  mass  Mg  in  time  dt  is  equal  to  the  increase  in  the  kinetic 
energy  of  the  mass  dMg  encompassed  by  the  shell  in  the  time  dt.  For  the 
case  of  an  extremely  porous  medium,  this  relation  becomes 

d_ 

‘  dt 

which  can  be  put  in  the  form 


Integration  indicates 


MR  *2  dM 
(_s _ \  R _ s 

(  2  >  "  2  dt 


d(MgR) 

dt 


=  0 


(35) 


M  R  =  CONST 
s 

M  R2  ,2 

s  _  (CONST)  ^ 

2  2M  ~  Y 

s 


(36a) 


(36b) 


The  constant  of  integration  is  undefined  in  Ref.  2.  Equations  (36a)  and 
(36b)  indicate  conservation  of  shell  momentum  and  a  decrease  in  shell 
kinetic  energy  with  increase  in  Y.  The  present  analysis,  Eq.  (33),  indi¬ 
cates  an  increase  in  momentum  and  a  slower  decrease  of  kinetic  energy  with 
increase  in  Y.  The  difference  is  due  to  the  fact  that  the  present  analysis 
includes  the  effect  of  gas  pressure  on  shell  momentum  [Eq.  (26)),  while 
Ref.  2  does  not,  [Eq.  (35)].  This  difference  is  a  consequence  of  the  pre¬ 
sent  assumption  that  the  shell  bounds  the  gas  during  the  entire  expansion. 
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III.  CONCLUDING  REMARKS 


The  expansion  into  vacuum  of  a  gas  bound  by  a  shell-like  mass  with 
zero  tension  has  been  evaluated.  The  flow  is  self-similar J  The  effect  of 
the  shell  is  to  reduce  expansion  velocities,  reduce  gas  dynamic  pressure, 
and  increase  the  net  impulse  at  a  station. 

An  approximate  local  similarity  solution  has  been  presented  for  the 
case  where  the  shell  is  surrounded  by  a  highly  compressible  porous  medium. 
The  shell  mass  at  each  instant  during  the  expansion  is  assumed  to  equal  the 
original  shell  mass  plus  the  mass  of  the  displayed  porous  medium.  In  the 
limit  R/Rq  >>1,  the  net  shell  momentum  increases  with  increase  in  R/Rq. 

The  latter  result  differs  from  the  corresponding  solution  for  a  point 
explosion  in  a  porous  medium  in  Ref.  2  due  to  inclusion  herein  of  the 
effect  of  gas  pressure  on  late  time  shell  motion. 
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APPENDIX  A 
SYMBOLS 


aQ  speed  of  sound  at  (r,t)  =  (0,  0) 

C  constant  related  to  M./M,.,  Eq.  (8c) 

6  s 

Eg  initial  energy  of  gas 

El  normalized  energy  integral,  Eq.  (8b) 

Iq  net  impulse  at  plane  of  synnetry  at  r  =  0,  Eq.  (14) 

Ig  impulse  per  unit  area  at  r1  >>  rQ  due  to  gas  flow,  Eq.  (21) 

Ig  impulse  per  unit  area  at  r^  >>  Tq  due  to  shell,  Eq.  (22) 

Jt  Ig  +  Js 

Mg  mass  of  gas 

Mg  mass  of  shell 

MI  normalized  mass  integral,  Eq.  (8a) 
p  gas  pressure 

ps  gas  pressure  at  n  =  1 

p0  gas  pressure  at  (r,t)  =  (0,0) 

R  shell  location,  R(t) 

Rq  initial  shell  location 

r  radial  distance 

rQ  intial  edge  of  gas  cloud,  rQ  =  Rq 
radial  distance  satisfying  r^  >>  Tq 
t  time 

v  velocity  in  r  direction 

Y  ratio  specific  heats 

6  2,  2tt  and  4*  for  a  =  0,1  and  2 

n  similarity  variable,  r/R 

p  gas  density 

pQ  gas  density  at  (r,t)  =  (0,0) 
pp  density  of  porous  medium 

a  0,1  and  2  for  planar,  cylindrical  and  spherical  flows 
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Subscript 


0  conditions  at  (r,t)  =  (0,0)  or  at  t  =  0 
1  condition  at  r1  >>  rQ 

Superscript 

(')  d(  )/dt 
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APPENDIX  B 

INTEGRALS  FOR  MI  AND  El 


Integral  expressions  for  MI  and  El  are  given  by  Eqs.  (8a)  and  (8b). 
These  have  been  integrated  in  closed  form.  The  results  for  o  =  1  are  given 
in  Eqs.  (9b)  and  (10b).  The  results  for  o  =  0,2  are  given  herein. 

The  integral  expression  for  MI  and  El  can  be  put  in  the  form 


1  ?  (N-1)/2 

MI;  El  =  (o+l)  J  (1-Cn  )  n°dn  (B-la) 

0 

where,  for  MI 


N  =  (y+1)/(y-1) 


(B-lb) 


and  for  El 


N  =  (3y-1)/(y-1) 


(B- 1c) 


For  y  =  5/3,  7/5,  9/7,...  the  parameter  N  has  the  values  4,6,8,...  and 

6,8,10,...  for  MI  and  El,  respectively.  The 
1/2 

substitution  C  n  =  sine  yields 


MI;  El  = 


o+ 1 


sin-V/2 


c(o*1)/2  J 


(cos  e 


o  N+2.  . 

2  cos  )d0 


(B-2) 


Since  N  is  an  even  integer  for  cases  of  interest,  Eq.  (B-2)  can  be 
integrated  to  yield  (for  o  =  0,2) 
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o+ 1  ( .  a  N+1 \r ( 1-C) 

MI;  EI  =  ~ 2  l1  "  5  uI5il - 


(N-1)/2 


2  N+2' 


N 


N-1  (1-C)(N"3)/2 
N  N-2 


N-1  m-3  (1-C)(N~5)/2  N-1  N-3  •••  5  (1-C)3/2 

N  N-2  N-4  +  *•'  +  N  N-2  •••  6  4 


N-1  N-3  N-5  •••  3  t( 1-C)1/2  sin"1  C1/2i 
N  N-2  N-4  •••  4  1  2  +  2  1/2  J 


0+1  o  ( 1-C) 


(N+n/2 


,0/2  2  N  +  2 


The  corresponding  value  for  o  =  1  is 


«:  EI  *  TBTTTc  I1  -  ('-«<N*,)/2I 


(B-3a) 


(B-3b) 


In  the  limit  C  =  1 


MI;  EI 


_  _2_ 
=  N+1 


(o  =  1) 


( B-4a ) 


=  (o+1)[l 


o  N+1 1 [N-1  N-3  •••  3 1 
2  N+2Jl  N  N-2  • ••  4j  4 


(o  =  0,2)  ( B-4b ) 
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